Two-Phase Regime in the Magnetic Field-Temperature Phase Diagram of a Type-II 

Superconductor 
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The magnetic field and temperature dependencies of the magnetic moments of superconducting 
crystals of VaSi have been studied. In a constant magnetic field and at temperatures somewhat below 
the superconducting transition temperature, the moments are hysteretic in temperature. However 
the magnetic moment-magnetic field isotherms are reversible and exhibit features that formally 
resemble the pressure-volume isotherms of the liquid-gas transition. This suggests the existence of 
a first-order phase transition, a two-phase regime, and a critical point in the superconducting phase 
diagram. The entropy change, determined from the data using the Clausius-Clapeyron equation, is 
consistent with estimates based on the difference in the vortex densities of the two phases. 
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The study of magnetic field-temperature phase dia- 
grams of type-H superconductors has revealed a rich 
range of behaviors, with the various transitions between 
vortex configurations, i.e., from lattice to liquid, glass to 
liquid, etc., resembling those found in ordinary matter. 
As a consequence, the study of thermodynamic proper- 
ties of "vortex matter" has emerged as an important field 
of research [EJ . The richness of the physics became appar- 
ent only after high-quality single-crystal samples of high 
temperature superconductors became available Q . Con- 
ventional superconductors that were subsequently stud- 
ied, were also found to exhibit related phenomena, ex- 
amples of which are phase transitions Q associated with 
the peak effect [Q in Nb, and coexistent vortex phases 
in the peak effect regime of NbSe 2 j§. The work on 
Nb combined magnetic measurements with small angle 
neutron scattering studies, and that on NbSe2 involved 
a direct map of the local ac induction using a variant 
of scanning Hall-probe microscopy. Domains of differing 
pinning strength remained distinct during the transfor- 
mation from stronger pinning and higher critical current 
at high temperature to weaker pinning and lower critical 
current at low temperature. Such behavior is a hallmark 
of a first order transition. In this Letter we report mea- 
surements of the equilibrium magnetic moments of VaSi 
crystals as a function of magnetic field and temperature. 
We provide an interpretation of the results using an anal- 
ogy with the phase diagram of the liquid-gas system that 
further supports the existence of coexistent phases in ad- 
dition to the existence of a critical point. Our approach, 
being based on a thermodynamic analysis of magnetiza- 
tion data alone, does not reveal the nature of the phases. 

The intermetallic compound VsSi is a type-II super- 
conductor that has been known since the 1950s §]. Sto- 
ichiometric and near-stoichiometric samples form in the 
A15 crystal structure. The crystals used in the present 
work were prepared using an electron-beam floating zone 



technique . They were characterized by metallographic 
and lattice parameter measurements. Also their resistiv- 
ity ratios, superconducting transition temperatures, and 
structural transition temperatures were cataloged. The 
samples were flat wafers cut from the VsSi rods pro- 
duced in the floating zone process, and polished using 
standard abrasives. Superconducting transition temper- 
atures ranged from 16.52 K to 16.97 K, and resistivity ra- 
tios from 2.6, for the lowest transition temperature sam- 
ple, up to 30 to 40 for samples with transition temper- 
atures the order of 16.95 K. Lattice parameters ranged 
from 4.7300 A for disordered samples down to 4.7242 A 
for the most ordered ones. Although some wafers were 
single crystals, most of them were polycrystalline with 
large crystallites. 

Measurements of magnetization were made using a 
Quantum Design Magnetic Properties Measurement Sys- 
tem (MPMS). Samples, with their flat sections parallel to 
the magnetic field direction, were attached to a stepper- 
motor controlled platform and driven through a super- 
conducting detection coil, which is a second-order gra- 
diometer connected to a superconducting quantum inter- 
ference device (SQUID). Standard software was used to 
determine the magnetic moment from the SQUID out- 
put signal. Although this is a conventional system for 
the determination of the properties of superconducting 
and magnetic materials, some detail is relevant, as the 
on-axis field homogeneity of the particular system (serial 
#1 of the MPMS instruments) is only 4 parts in 10,000 
over a distance of 3 cm. This means that during a 3 cm 
scan at a field of 15 kOe, the sample in effect will be sub- 
ject to a variation of magnetic field of order 6 Oe. Much 
has been written about the effect on the reversibility of 
measurements of field inhomogeneity in this instrument 
||[)| . We do not believe that much of this is relevant for 
this particular set of studies with these samples as the 
m vs. H curves that were obtained were reversible under 
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FIG. 1. Three dimensional plot showing the magnetic mo- 
ment measured at fixed temperature, as a function of applied 
magnetic field and temperature. The temperature ranges 
from 13.8 K to 15.2 K, in 0.1 K increments. 

most conditions. It is likely that the variation of the field 
during a scan is too small a fraction of the total field to 
introduce any significant additional irreversibility. 

Standard plots of the magnetic moment vs. magnetic 
field, m(H), for these VaSi wafers were consistent with 
the usual results for type II superconductors, with H c \ 
measured to be 700 Oe. Because the upper limit on the 
magnetic field in this apparatus was 55 kOe, and the 
upper critical field of VsSi at low temperatures is 200 
kOe, there was a relatively high lower bound on the tem- 
perature at which a complete set of data could be ob- 
tained. In Fig. |l| we show a three-dimensional plot of 
raw data of m(H) at different fixed temperatures, for a 
temperature regime relatively close to the zero-field criti- 
cal temperature. This particular sample was a wafer that 
was 0.76 cm x 0.48 cm in area and 0.053 cm thick. Be- 
low a characteristic temperature, flat sections appear in 
m(H) isotherms, and as temperature is decreased further 
these deform into structures which appear to oscillate re- 
versibly (without hysteresis). The flat and oscillatory 
regimes separate regions which are distinguished by dif- 
ferent values of the susceptibility, (dm/dH) T . In Fig. || 
we show m(T) at fixed H , for a range of magnetic fields 
spanning the regime of T and H in which there are non- 
trivial features in the m(H) isotherms. For this mode 
of data acquisition, hysteretic behavior in temperature 
is clearly observed. This occurs even though there is no 
apparent hysteresis in the m(H) isotherms over the same 
range of temperatures. 

The above-described behavior was observed in five out 
of seven VaSi samples. It was always seen if the samples 
are cooled in zero field to low temperatures, and the field 
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FIG. 2. Magnetic moment as a function of temperature at 
different applied fields: 25 (bottom), 30, 35, 40, 45, 50, and 55 
x 10 3 Oe (top). Closed symbols represent data taken as the 
temperature is increased, while open symbols represent de- 
creasing temperature data. Noticeable hysteresis is observed 
as the field is increased beyond 35 x 10 3 Oe. 

and m(H) curves are obtained with increasing tempera- 
ture. On the other hand, reversibility was lost in most 
of the samples if they were cooled in high fields before 
a measurement cycle was initiated. Even in this circum- 
stance curves of m(H) in increasing field were the same 
as those obtained on cooling in zero field. 

It is apparent from the m(H) isotherms shown in Fig. |l| 
that the thermodynamics of the system may best be de- 
scribed by considering the magnetic work dW done by the 
system when H is varied at constant m, dW — —mdH 
|l0| . This also corresponds to the choice in which the 
internal energy does not include the field energy. The 
magnetic moment, m, in this instance is analogous to the 
pressure P whereas the magnetic field, H, is analogous 
to the volume V. The variables P and V are used to de- 
scribe mechanical work in liquids. Thus m(H) isotherms 
are analogs of P(V) isotherms. The data exhibit features 
that resemble those found in liquid-gas transitions. There 
are isotherms in which there are regimes in which m is 
independent of H . As the temperature is lowered, the 
"flat" region is replaced by a regime in which m rises and 
falls. This is reminiscent of the behavior of the van der 
Waals equation for a fluid, which can describe metastable 
persistence of the low and high density regimes, into 
what is usually a two-phase region, when the pressure 
is constant over a range of volumes. Only in this in- 
stance the magnetic moment appears to move smoothly 
and reversibly between the phases, without abrupt nucle- 
ation phenomena. Perhaps the appearance, with decreas- 
ing temperature, of a continuous transformation between 
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FIG. 3. Phase diagram in H — T space. The upper symbols 
denote H C 2(T). The shaded area denotes the putative two 
phase region bounded by the coexistence curve (diamonds). It 
is not known whether this region extends to low temperatures 
because of the limitation of magnetic field in the measuring 
instrument. 



the two states, has to do with the growth of barriers to 
changing the vortex configuration and moving vortices. 

The existence of flat isotherms and the analogy with 
van der Waals isotherms for the oscillatory isotherms sug- 
gests a two-phase regime with an endpoint H]. The 
analogy with the fluid system suggests that there is a first 
order transition line terminating at a critical point. The 
two phases associated with the first order phase tran- 
sition would have to have the same symmetry |T^ ]. In 
Fig. U we show a phase diagram for the behavior of this 
sample. The figure includes H C 2(T), in addition to lines 
denoting the high and low-field limits of the two-phase 
regime, delineated by the shaded region. Because of the 
limitation to 55 kOe on the highest field achievable with 
this apparatus, we were not able to explore this behavior 
to lower temperatures than shown in Fig. [l]. As a con- 
sequence we do not know whether the two-phase regime 
terminates at a second critical point as the temperature 
is reduced, or persists to arbitrarily low temperatures. 

The results discussed here are very different from those 
that have been reported previously for VaSi 0,0. 
Their m(H) isotherms were always hysteretic, exhibit- 
ing peaks and dips with decreasing and increasing mag- 
netic field, respectively. Similar behavior was found in 
our most disordered samples, in our thick samples, and 
in our less disordered samples when their flat parts were 
rotated such that they were at an angle with respect to 
the field direction. The details of m(H) curves would 
thus appear to depend critically on the geometry of mea- 



surement, and the degree of disorder, although a detailed 
mechanism is not known. (The fact that the samples in 
which these effects were seen were very thin wafers may 
account for our results as this might minimize demag- 
netizing effects as well as the role of off-axis field inho- 
mogeneities.) We must emphasize that our samples also 
have dislocations and defects, and are usually not single 
crystals, so our observations cannot be a consequence of 
having samples of "superior" quality. We believe that 
some combination of defects and geometry allow a two- 
phase equilibrium state to be achieved for some, but not 
for the full range of parameters. 

There have been many discussions of nonlinear features 
of type-II superconductors. In fact there has been some- 
what of a debate as to whether in paramagnetic materials 
such as VaSi these features are a signature of a general- 
ized Fulde-Ferrell- Larkin-Ovchinnikov (GFFLO) phase 
|l5f which can appear at high rather than at low tem- 
peratures as originally predicted. In this phase, the su- 
perconductivity becomes spatially inhomogeneous, with 
a wavelength that diverges at the transition to the normal 
state. The transition to the GFFLO phase with increas- 
ing temperature is predicted to be first order. There are 
no obvious features of the data that would support this 
interpretation. 

An alternative explanation is that our observations are 
associated with the peak effect, a phenomenon known 
for many decades in conventional superconductors [0. 
In the peak effect, the critical current, instead of de- 
creasing as the superconductor-normal phase boundary is 
approached, increases before falling to zero at the tran- 
sition. The reversibility of our m(H) isotherms would 
seem to rule out our observations as being attributable 
to dynamical effects associated with the critical current. 
Our previous discussion points to two distinct equilib- 
rium vortex phases, possibly glass phases, or disordered 
lattice phases separated by a first order transition, which 
terminates at a critical point. It is conceivable that in a 
dynamical regime these phases could have different pin- 
ning currents as reported in Ref . p) . 

An important issue is the question of what distin- 
guishes the two phases whose thermodynamic signature 
is inferred in our interpretation of m{H) isotherms. The 
high- and low-field regimes on either side of the "two- 
phase" regime have different magnetic susceptibilities 
(dm/ dH) T and their vortex densities are different. Thus 
the two coexisting phases would have different degrees of 
screening as m is fixed with H changing, so that B would 
have to change. In equilibrium the free energies of the two 
phases at constant T and M would be the same. With 
our choice of dependent and independent variables, this 
leads to a Clausius-Clapeyron equation of the form 
dm/dT = AS /AH. 

A thermodynamic consistency check of a picture in 
which the two phases are distinguished by their vortex 
densities is possible. One can determine AS from mea- 
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sured quantities and compare this value with an estimate 
obtained from the change in the number of states asso- 
ciated with the change in the vortex density across the 
transition. To obtain AS from the isotherms, we first 
estimate dm/dT from the data of Fig. |l|. (For the oscil- 
latory isotherms we use the average moment.) We find 
dm/dT to be 0.005 emu/K, and relatively constant over 
the accessible range of temperatures. The change in field 
across the "flat" region is about 3000 Oe. Using these 
values A5 is found to be approximately 10 17 fcs. 

Then, the entropy per configuration in the vortex state, 
S/ks, is given by ln[iV e ], where the number of configu- 
rations N e w N V N C . Here N v is the number of vortices 
threading the sample, and N c the number of configura- 
tions of a single vortex. The number of vortices N v is 
given by the flux divided by the flux quantum. Taking 
the relevant area of the crystal to be 0.05 x 0.5 cm 2 , 
and the average magnetic field to be 40,000 Oe, we find 
N v = 5 x 10 9 vortices. The number N c can be esti- 
mated as the number of vibrational modes of a vortex 
line. Thus iV c is of the order of the number of unit cell 
layers, which with a sample of length 0.7 cm, is 1.4 x 10 7 
unit cell layers. Then the change in entropy associated 
with the change in field between values 1 and 2 is given 
by S 2 - Si « k B [N v2 - N vl ] N c In (N av N c ) . Here N v2 
and N v i are the numbers of vortices associated with the 
two fields, and N av is the average number of vortices in 
the two-phase regime. With an average field of 40,000 
Oe, and a field change of 3,000 Oe in the case of the 
data in question, S 2 — Si is found to be « 2.2 x 10 17 ks 
in remarkable agreement with the value estimated from 
the data. If the coherence length, £ (0) , estimated from 
dH c2 /dT to be approximately 27 'A, determines the cutoff 
in the number of modes, this estimate would be a factor 
of five lower. Estimates of the entropy change taking 
into account only the entropy of the vortex cores yield a 
somewhat smaller estimate of the entropy change. 

In summary, magnetic moment -magnetic field 
isotherms of crystals of superconducting VaSi have been 
found to resemble the pressure-volume isotherms of the 
liquid-gas transition. We have suggested that this is ev- 
idence of a first-order phase transition and a two-phase 
regime, with a critical point in the superconducting phase 
diagram. As mentioned, in this instance, the two phases 
would have to possess the same symmetry. The pre- 
cise nature of vortex phases in the presence of weak dis- 
order is quite complicated |Tj]|, so that further studies 
such as measurements of specific heat, ac magnetic sus- 
ceptibility, and small angle neutron scattering would be 
needed to resolve this issue. This work was supported 
in part by the National Science Foundation under grant 
NSF/DMR- 0138209. 
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